Accurate numerical solution to the finite-size Dicke model 
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By using extended bosonic coherent states, a new technique to solve the Dicke model exactly 
is proposed in the numerical sense. The accessible system size is two orders of magnitude higher 
than that reported in literature. Finite-size scaling for several observables, such as the ground- 
state energy, Berry phase, and concurrence are analyzed. The existing discrepancy for the scaling 
exponent of the concurrence is reconciled. 



PACS numbers: 03.65.Ud,03.67.Mn,42.50.-p,64.70.Tg 

The Dicke model describes the interaction of N two- 
level atoms (qubits) with a single bosonic mode and has 
been paradigmatic example of collective quantum behav- 
ior. It exhibits a " superrandiant" quantum phase tran- 
sition (QPT)[2| in the thermodynamic limit, which was 
first studied by Hepp and Lieb Q at weak coupling. In 
recent years, the Dicke model has attracted consider- 
able attentions due to the progress in the QPT. First, 
entanglement, one of the most striking consequences of 
quantum correlation in many-body systems, shows a deep 
relation with the QPT0|. Understanding the entangle- 
ment is also a central goal of quantum information sci- 
ence. Second, Berry phase has been recently drawn in 
the study of the QPT A drastic change at critical 
point in the QPT may be reflected in the geometry of 
the Hilbert space, the geometric phase may capture the 
singularity, and therefore signal the presence of the QPT. 

In the thermodynamic limit, the Dicke model is ex- 
actly soluble in the whole coupling range, based upon 
the Holstein-Primakoff transformation of the angular mo- 
mentum algebra 6]. For finite N, the Dicke model is in 
general nonintegrable. The finite-size correction in this 
system has been shown to be crucial in the understanding 
the properties of the entanglement @, H S E3, El E3 and 
the Berry phase 13] , in which one can characterize univer- 
sality around the critical point in the QPT. However, a 
convincing exact treatment of the finite-size Dicke model 
is still lacking. To the best of our knowledge, the finite- 
size studies are limited to numerical diagonalization in 
Bosonic Fock state t 6][7][9] in small size system N < 35, 
the adiabatic approximation and l/N expansion 

based on modified Holstein-Primakoff approach [111]. 

Recently, the Dicke model is closely related to many 
fields in quantum optics and condensed matter physics, 
such as the superradiant behavior by an ensemble of 
quantum dots [3] and Bose- Einstein condensates [l5|, 
and coupled arrays of optical cavities used to sim- 
ulate and study the behavior of strongly correlated 
systems [l^]. The finite-size Dicke model itself is also of 
practical interest. It was observed that the Dicke model 
for finite N can be realized in several solid-state systems. 



One Josephson charge qubit coupling to an electromag- 
netic resonator [TJj can be described by N — 1 Dicke 
model , which is just special case of spin-1/2 (Rabi Hamil- 
tonian). The features of the superconducting quantum 
interference device coupled with a nanomechanical res- 
onator may be captured in the framework of the finite- 
size Dicke model (la ]. 

In this paper, we propose an exact technique to solve 
Dicke model numerically for finite N by means of ex- 
tended bosonic coherent states. The correlations among 
bosons are added step by step until further corrections 
will not change the results. Some important quantities 
are calculated exactly. 

Without the rotating-wave approximation, the Hamil- 
tonian of N identical qubits interacting with a single 
bosonic mode is originally given by 



2A 

H = cua + a + AJ Z H — i={ a+ + a )Jx, 



N 



(1) 



where a + and a are the field annihilation and creation 
operators, A and u> are the transition frequency of the 
qubit and the frequency of the single bosonic mode, A is 
the coupling constant. J x and J z are the usual angular 
momentum. There is a conserved parity operator n = 



a m(Jz+N/2+a+a) 



which commutes with the Hamiltonian 



(1). For convenience, we use a transformed Hamiltonian 
with a rotation around an y axis by an angle ? , 

A 2A 

H = uja + a-— (J+ + J-) + —={a + + a) J z , , (2) 



N 

where J+ and J_ are the angular raising and lowing 
operators and obey the SU(2) Lie algebra [J+,J_] 
2.1 [J Z ,J±\ = ±J±. So the Hilbert space of this alge- 
bra is spanned by the Dicke state { \j, m) , m = —j, —j + 
1, ...j — 1, j} with j — N/2, which is the eigenstate of J 2 
and J z with the eigenvalues j(j + 1) and m. 

The Hilbert space of the total system can be expressed 
in terms of the basis {\<fm} b ® \j , m}} , where only the 
state of bosons|v? m )t, is to be determined. A "natural" 
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basis for bosons is Fock state \l) b — [(a + )*/v/!] |0) h . In 
the Dicke model, the bosonic number is not conserved, so 
the bosonic Fock space has infinite dimensions, the stan- 
dard diagonalization procedure (see, for example, Ref. 
[@00]) i s t° apply a truncation procedure considering 
only a truncated number of bosons. Typically, the con- 
vergence is assumed to be achieved if the ground-state 
energy is determined within a very small relative errors. 
Within this method, one has to diagonalize very large, 
sparse Hamiltonian in strong coupling regime and/or in 
adiabatic regime. Furthermore, the calculation becomes 
prohibitive for larger system size since the convergence of 
the ground-state energy is very slow. Interestingly, this 
problem can be circumvented in the following procedure. 

By the displacement transformation A m — a + g m with 
g m = 2\m/oj^/~N \ the Schrb dinger equation can be de- 
scribed in columnar matrix, and its m row reads 

- A Jm \<Pm) b \j, m + 1) - Aj~ \<p m ) b \j, m-1) 
uj(A^A m - g 2 m ) \<p m ) b \j, m)=E \<p m ) b \j, m) , (3) 

where = \\J{j{j + 1) — m(m ± 1). Left multiplying 
{(n, j\\ gives a set of equations 

-Aj+ Wn+i) b - Aj- \<p n -i) b + u{A+A n - g 2 n ) \<p n ) b 
= Etp n \<p n ) b , (4) 

where n = —j, —j + 1, ...j. 

Note that the linear term for the bosonic operator 
a(a + ) is removed, and a new free bosonic field with op- 
erator A(A + ) appears. In the next step, we naturally 
choose the basis in terms of this new operator, instead of 
a(a + ), by which the bosonic state can be expanded as 

N tr 

Wn) b = £ C ^(A+) fe |0) A „ 
fc=0 



^2 C n,k-7=(a + +g n ) 
k=0 



|0) B ,(5) 



where N tr is the truncated bosonic number in the Fock 
space of A(A + ). As we know that the vacuum state |0)^ 
is just a bosonic coherent-state in a(a + ) with an eigen- 
value On [T9L I20I |2"H . So this new basis is overcomplete, 
and actually does not involve any truncation in the Fock 
space of a(a + ), which highlights the present approach. It 
is also clear that many-body correlations for bosons are 
essentially included in extended coherent states (5). Left 
multiplying state An (l\ yields 



w(Z 



)Cnl- Aj r 



N tr 

E 

k=0 



Cn+x,k A n {l \k) J 



- Aj n ^2c n ^ 1:k An {l\k) An i ) = Ec n j, (6) 



where 



k=0 

A n (l\k) An _ 



(-l)'A.fc, 
(-l) k D Lk , 
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FIG. 1: The scaled ground-state energy Eo / (j A) as a function 
of A by DCS and DFS with different truncated numbers. The 
inset show the scaled ground-state energy as a function of the 
N for fixed A = 1. 



with 



min[l,k] , , 
= , G 2 /2 — - 
r=0 



(I - r)\(k - r)lr\ 



2A 



N 



(7) 



. Eq. (6) is just a eigenvalue problem, which can be 
solved by the exact Lanczos diagonalization approach in 
dimensions (N + l)(N tr + l)[22j|. The integral of motion 
due to the parity II can simplify the computation fur- 
ther. To obtain the true exact results, in principle, the 
truncated number Ntr should be taken to infinity. For- 
tunately, it is not necessary. It is found that finite terms 
in state (5) are sufficient to give very accurate results 
with a relative errors less than 10~ 6 , in the whole cou- 
pling range. We believe that we have exactly solved this 
model numerically. 

To show the effectiveness of the present approach, we 
first calculate the ground-state energy Eq by diagonal- 
ization in coherent-states (DCS) as described above , 
and compare with those by numerical diagonalization in 
Fock state (DFS) 0] for N=32. For convenience, we 
introduce two dimensionless parameters D = A/uj and 
a = 4A 2 /Aw. Fig. 1 shows the scaled ground-state en- 
ergy E /(jA) as a function of A when the Hamiltonain 
is on a scaled resonance A = oj = 1 . For the same trun- 
cated number N tr — 6, the ground-state energy by DCS 
is much lower than that by DFS. With the increase of 
Ntr, DFS results approach the DCS one monotonously 
in the stronger coupling regime. One can find that the 
ground-state energy by DCS with N tr = 6 is even much 
better than that by DFS with Ntr = 100. In addition, 
the DFS results become worse and worse as the N in- 
creases for fixed N tr and A, as clearly shown in the inset 
of Fig. 1. 

The ground-state energy is non-analytical at the crit- 
ical point of the QPT in the thermodynamic limit. To 
show the effect of the precursor of the QPT on the present 
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FIG. 2: The relative ground-state energy difference as a func- 
tion of -/V for different truncated numbers Ntr at critical point. 
The inset show 8E(N tr ) as a function of A at fixed N for dif- 
ferent Ntr- 



FIG. 4: Scaling of the ground-state Berry phase 1+ < J x > 
/N as a function of N at the critical point for D = 0.1, 1, 
and 5. The upper inset show the corresponding slope as a 
function 1/N, and the lower one gives the Berry phase versus 
a for different N at D = 5. 




iog 2 (N) 



FIG. 3: Scaling of the ground-state energy Eo/ND + 0.5 as 
a function of N at the critical point for D = 0.1, 1, and 10. 
The inset show the corresponding slope versus 1/N. 

approach, we display the relative ground-state energy dif- 
ference 6E(N tr )/E% = (E^ tr - E°)/E° as function of N 
at the critical point A c in Fig. 2. It is interesting to 
find that as N increases, the exact results can be ob- 
tained within smaller truncated number. This feature 
facilitates the calculation for rather large system size. In 
the inset of the Fig. 2, we can find that in both the weak 
and strong coupling regime, the results converge rapidly 
with the truncated number N tr for fixed large N, indict- 
ing that the Dicke model can be solved easily in these 
two limits. More interestingly, just in the critical regime, 
a dip structure is found in all curves, indicating a trace 
of QPT at finite size system. 

To study the universality of the QPT, we will compute 
the finite-size scaling exponents for several observables. 

Ground-state energy. -We calculate the ground-state 
energy for different values of D. In the thermodynam- 



ical limit, Eq/ND = —0.5. We are able to study the 
system up to N — 2 12 atoms in a PC. To show the lead- 
ing finite-size corrections, we plot Eq/ND + 0.5 versus 
N for D = 0.1, 1, 10 in a log-log scale in Fig. 3. The 
inset presents the slope as a function 1/N. It is observed 
that the slope of all these curves in the large N regime 
give a universal exponent —1.0 ± 0.02, from the above 
and below respectively, consistent with that by a mod- 
ified Holstein-Primakoff approach [Tlj]. The asymptotic 
behavior demonstrates that the exponent by no means 
lies on the vicinity of —4/3 in the the adiabatic regime 
reported in Ref. (l0| . 

Berry phase.- To generate the Berry phase, we intro- 
duce a time-dependent unitary transformation U(T) = 
exp[— icj)(t) J x ], where </>(t) is changed adiabatically and 
slowly from to 2ir. The ground-state Berry phase can 
be defined asp^j 

r 2rr d 
7 = i / < ip\U + (t) — U(t)\ip >d<j) = 2n< <p\J x \ip > 
Jo dp 

where \tp > is just the time-independent ground-state 
wavefunction, which can be obtained by the present DCS. 
In the thermodynamical limit < J x > /N = — 1. 

Fig. 4 shows the scaled ground-state Berry phase 1+ < 
J x > /N as a function of iV for different values of D in 
log-log scale. A power law behavior exists in the large 
N . One can see from the upper inset that the finite-size 
exponents extracted from all curves tend to a converging 
value (—0.67 ± 0.02). This result is in good agreement 
with that in Ref. ([l3j) with a adiabatic approximation. 
We find that the scaling exponent for the Berry phase in 
Dicke model is also universal, not merely in the adiabatic 
limit. The Berry phase tends to be discontinuous in the 
critical regime with the increase of N, shown in the lower 
inset. 
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FIG. 5: Scaling of concurrence as a function of N at the criti- 
cal point for D = 0.1, 1, and 5. The upper inset show the slope 
versus 1/N, the lower inset displays the scaled concurrence as 
a function of A/A c . 

Pairwise entanglement. -An important ground-state 
quantum essential feature is the concurrence, which 
quantifies the entanglement between two atoms in the 
atomic ensemble after tracing out over bosons. The 
scaled concurrence (entanglement) in the ground-state 
can be defined as C N = 1 - 4 < J% > /N [ill [1| . In the 
thermodynamic limit N — > oo, the scaled concurrence 
at critical point Coo(A c ) can be determined in terms of 
D Q. For comparison, we can calculate the quantity 
Coo (A c ) — Cat (A c ). In Fig. 5, we present this quantity as 
a function of N for different values of D in log-log scale. 
Derivative of these curves is presented in the upper in- 
set, and the exponent of concurrence is estimated to be 
—0.33 ± 0.02. For large D, the exponent for concurrence 



is only given reasonably in the asymptotic regime. From 
the lower inset, we find that a cusp is formed with in- 
creasing N in the critical regime, demonstrating the relic 
of the QPT in finite size. 

Lambert et al., performed a numerical DFS and give 
the exponent of concurrence —0.25 ± 0.01 [9j with N < 
32. Reslcn et al., derived a effective Hamiltonian, and 
found the exponent to be -0.26 ± 0.02 with N < 35. [13] 
Recently, Vidal et al. using the diagonalizing a expanded 
Hamiltonian at order 1/N based on Holstein-Primakoff 
representation, and predicted finite-size scaling exponent 
for concurrence 1/3 [ll|, agreeing well with our results. 
Due to the huge system size we can touch here, the pre- 
vious controversy can be attributed to the small system 
size @,[T3]. 

In summary, we have proposed a numerically exact so- 
lution to the Dicke model in the whole coupling regime 
for a huge system size up to N = 2000 ~ 4000, which is 
almost two orders of magnitude higher than that reported 
in literature [H, H3] ■ Several quantities related to QPT 
such as the ground-state energy, the Berry phase, and 
the concurrence are then calculated exactly. It is found 
that the scaling exponent for the ground-state energy is 
different from that in a previous study 0. The previous 
discrepancy for the exponent of concurrence is reconciled 
by the present calculation in very large size systems. The 
precise estimate of the scaling exponent for these quanti- 
ties is very significant to clarify the universality of QPT. 
The methodology sketched here is helpful to study a large 
number of spin(fcrmion)-boson coupling systems. 
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